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Abstract 

In many Lagrangian field theories one has a Poisson bracket de- 
fined on the space of local functionals. We find necessary and sufficient 
conditions for a transformation on the space of local functionals to be 
canonical in three different cases. These three cases depend on the 
specific dimensions of the vector bundle of the theory and the associ- 
ated Hamiltonian differential operator. We also show how a canonical 
transformation transforms a Hamiltonian evolutionary system and its 
conservation laws. Finally we illustrate these ideas with three exam- 
ples. 
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1 Introduction 

In Hamiltonian evolutionary systems one has a Poisson bracket. This Poisson 
bracket is defined on the space of local functionals of a fiber bundle. The 
system is characterized by a Hamiltonian which is a local functional and the 
associated Poisson bracket. The Poisson bracket in turn is characterized by 
a Hamiltonian differential operator (see the discussion in the next section). 
The solutions to the evolutionary system are sections of the fiber bundle. 
Generally an evolutionary system has the form 



ut = V5H 



where T> is the associated Hamiltonian differential operator, 5 the variational 
derivative and Ti the Hamiltonian of the system. The variational derivative 
of a local functional V = J M Pv has the same components as E(Pu) where 
E is the Euler-Lagrange operator and P is a local function representing the 
local functional V, while v is a volume element on the base manifold M. A 
good reference on evolutionary systems and their theory can be found in j3] 
and |B] for example. The reader may consult [U] as well. 

An action on the bundle induces an action/transformation on the space 
of local functionals. If this transformation preserves the Poisson bracket we 
call it a canonical transformation. This terminology was used in PQ and is 
in anology to the one used in [7j for the case of symplectic manifolds. In 
PP canonical transformations were studied in the case the Poisson bracket 
is defined by a differential operator of order zero. In this paper we study 
canonical transofrmations in the case of higher order differential operators. 
We work with vector bundles with m-dimensional fibers and n-dimensional 
base space. We consider special cases of higher order differential operators 
which we find applications for and intend to broaden our study in future work. 

A canonical transformation will transform the evolutionary system, its 
Hamiltonian, and its conservation laws where a new system is obtained. We 
will show that the conditions for a canonical transformation are rather re- 
strictive in some cases. We illustrate by a few examples at the end of the 
paper. 

It would be interesting to study the sh-Lie algebra structure and reduction 
as was done in jT] and [3]. However we leave this for a possible future work. 

2 Background material 

Let 7r : E — > M be a vector bundle of dimension m + n and with a base space 
an n-dimensional manifold M. Let J°°E be the infinite jet bundle of E. The 
restriction of the infinite jet bundle over an appropriate open set U C M is 
trivial with fiber an infinite dimensional vector space V°°. The bundle 

tt°° : J°°Eu = U xV°° 
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then has induced coordinates given by 

We use multi-index notation and the summation convention throughout the 
paper. If is the section of J°°E induced by a section <fi of the bundle E, 
then u a o = u a o ip and 

U }or<f ) = (d il d i2 --Ar)(u a or<j>) 

where r is the order of the symmetric multi-index / = i 2 , i r }, with the 
convention that, for r = 0, there are no derivatives. For more details see j2] 
and 0. 

Let Loce denote the algebra of local functions where a local function 
on J°°E is defined to be the pull-back of a smooth function on some finite 
jet bundle J P E via the projection from J°°E to J P E. Let Loc° E denote the 
subalgebra of Loce such that P G Loc° E iff (j°°0)*P has compact support for 
all <p G YE with compact support and where TE denotes the set of sections of 
the bundle E — > M. The de Rham complex of differential forms Q*(J°°E, d) 
on J°°E possesses a differential ideal, the ideal C of contact forms 9 which 
satisfy (j°°<fi)*9 = for all sections <fi with compact support. This ideal is 
generated by the contact one-forms, which in local coordinates assume the 
form 9 j = duj — ufjdx 1 . 

Now let Co denote the set of contact one-forms of order zero. Contact one- 
forms of order zero satisfy {j l <f))*(9) = and, in local coordinates, they as- 
sume the form 9 a = du a -u\dx\ Notice that both C and fT- 1 = ^{J^E) 
are modules over Loce- Let Qq' 1 denote the subspace of fi™' 1 which is lo- 
cally generated by the forms {(9 a A d n x)} over Loce- Let v denote a vol- 
ume element on M and notice that in local coordinates v takes the form 
v = fd n x = fdx 1 A dx 2 A ... A dx 11 for some function f : U R and U is a 
subset of M on which the af's are defined. 

d d 

Define the operator Di (total derivative) by Di = Tr7 + M a~Q — a ( reca U we 

C/X C/lij 

assume the summation convention, i.e., the sum is over all a and multi-index 
J). For I = ■• -ik} where k > 0, Dj is defined by Di = o D i2 o 
• • ■ o Di k . If / is empty I — {} then Di is just multiplication by 1. We also 
define (—D)j = (— 1)' 7 'D/. Recall that the Euler-Lagrange operator maps 
tt^iJ^E) into tt^iJ^E) and is defined by 

E(Pu)=E a (P)(9 a Au) 
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where P G Loce, v is a volume form on the base manifold M, and the 
components E a (P) are given by 



For simplicity of notation we may use E(P) for E(Pz/). We will also use Pj 
for — + wj-j— and E a (P) for (-P) 7 (_) so that E(P) = E a (P)(6 a A u) 
in the w a ) coordinate system. 



Now let JF be the space of functionals where V G T iff P = / P^ for 

some P G Loc° E . Let D be a differential operator which has components 
P ab = u abI Dj where uj abI G Loce and P/ is a combination of the total 
derivatives as determined by the multi- index J, i.e., Dj is a composite of the 
form Pjj o D i2 o ■ ■ ■ o D ik . Define a Poisson bracket on the space of local 
functionals T by 

{P,Q}(0) = / [E(P)V(E{Q))oj<f>]u, 

J M 



where G TP, v is a volume form on M, V = / Pi/, Q = Qv, and 

P,Q £ Loc° E . In this expression one may represent P as a matrix differential 
operator and E as a row/column vector (as appropriate) consisting of the 
components E a . We assume that P is Hamiltonian so that we have a genuine 
Poisson bracket that is antisymmetric and satisfies the Jacobi identity (e.g. 
see |8J. Using local coordinates {x^^uf) on J°°E, observe that for G TE 
such that the support of lies in the domain Q of some chart x of M, one 
has 

{P,Q}(0) = / ([E a (P)V ab (E b (Q))}o ^oj- 1 )^) 

where x~ l is the inverse of x = (x M ). 

The functions P and Q in our definition of the Poisson bracket (of lo- 
cal functionals) are representatives of P and Q respectively, since gener- 
ally these are not unique. In fact T ~ P™(J°°P), where H^(J°°E) = 
Q™> (J co E)/(imd H f)tt™'°(J oo E)) and imd H is the image of the differential 
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dji defined by du = dx % Di. We refer the reader to [T] for more details and 
for the notation. The interested reader may also consult and j3] for more 
on the de Rham complex and its cohomology. 

Let xp : E — > E be an automorphism, sending fibers to fibers, and let 
ipM '■ M — ► M be the induced diffeomorphism of M. Notice that ip induces 
an automorphism jip : J°°E — > J°°E where 

m{{r<t>){p)) = ° <p o ^-/)(^ M ( P )), 

for all cf) G YE and all p in the domain of 0. In these coordinates the indepen- 
dent variables transform via x M = ip^{x v ). Local coordinate representatives 
of i/jm and jip may be described in terms of charts (Q, x) and (Q, x) of M, 
and induced charts ((vr 00 )" 1 ^), (z**,^)) and ((vr 00 )- 1 ^), (#*,uf)) of J°°£. 

Observe that the total derivatives satisfy 

A(^o#) = {{DjF) ojip)Diip 3 M (2.1) 

where x- 7 = ^^(rr). 

Example Let M = R 2 , .E = R 2 x R 1 and consider the transformation ip 
defined by u = xu + 3yu 2 , x = x cos 9 + y sin 9,y — —x sin # + y cos Now let 
F = u x then F o jip = {xu x + u + 6yuu x ) cos 9 + [xu y + 3m 2 + 6yuu y ) sin 6* so 
that D X (F o j0) = [xu xx + 2m x + 6yu 2 + 6yuu xx ) cos + (xu yx + u y + 6uu x . + 
6yu x u y + 6yuu yx ) sin9. On the other hand the right-hand side of equation 
(|2.1|) yields [u xx cos9 + u xy (— sm9)} ojtp which when evaluated and simplified 
yields the same expression as above. 

Remark For simplicity we may skip writing the tilde's, so in the above 
example one simply writes F = u x instead of F = u x , ...etc. 

3 Canonical transformations of the Poisson 
structure 

Let L : J°°E — > R be a Lagrangian in Loce (generally we will assume 
that any element of Loce is a Lagrangian). Let L = L o (x m ,m") _1 and let 
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L = L o (x^Mj) 1 . Then, in local coordinates, L is related to L by the 
equation 

(LojYOdet(J)=L, 

where jip = [x v , w^) o (a;' 4 , and J is the Jacobian matrix of the transfor- 
mation ipM — x u o(x tl )~ 1 . With abuse of notation we may assume coordinates 
and charts are the same and write x v = iPm{%^)- For simplicity, we have also 
assumed that ipM is orientation-preserving. In this case the functional 

C= [ Ld n x 
Jq 

is the transformed form of the functional 

£= [ Ld n x 
Jq 

where L and L are related as above, Q is the domain of integration and Q 
is the transformed domain under jtp (see jHj pp. 249-250). Notice that both 

of these are local coordinate expressions of the equation £ = Lu, for 

J M 

appropriately restricted charts. Now suppose that ip is an automorphism of 
E, jijj its induced automorphism on J°°E, and ipM its induced (orientation- 
preserving) diffeomorphism on M. Also suppose that L and L are two La- 
grangians related by the equation (L o jip)det(ipM) = L. We have: 

Lemma 3.0.1 Let P be a Lagrangian as above, then 

E a ((Po^)det(^)) = det(Vivf)^f (E C (P) o (3.2) 

Proof First notice that E u a(L) = det(ip M )^^CE ii c(L)oj' l p) (see M pp.250). 

But (L o jip)det(ipM) = L. The identity 13.21 follows by letting P = L. 
Notice that this is justified since L is arbitrary in the sense that given any 
L' there exists an L derived from a Lagrangian L as above such that {JJ o 
ji/j)det(ipM) — L since jip is an automorphism. | 

Let ip denote the mapping representing the induced action of the auto- 
morphism on sections of E, i.e., t/> : YE — > TE where ip((f>) = ip o o and 
is a section of E. This induces a mapping on the space of local functionals 
given by 
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= / lP°j(^°<fr°^M)] U 
J M 

[Po ii)oj(\)o^)\v 

[P o jip o j(f)](detil>M)v, 
where 

and is a section of -E 1 . 

Now we find conditions on those automorphisms of the space of local 
functional under which the Poisson structure is preserved for a few special 
cases of the bundle E and differential operators V. 





3.1 Case I: dim(M) = l, dim(£)=2; T> = uD x 

We begin with the case where M is 1-dimensional, E is 2-dimensional and 
consider first order differential operators of the form V = uoD x with uo G 
Loce- For simplicity we will skip the tilde's in our notation. Recall that 

{V, Q}= [ P(E(Q))E(P)i/, and hence 
Jm 

{P,Q}(0)= / [V(E(Q))E(P)oj<j>]u. (3.3) 

In this case the Poisson bracket is preserved: {^(V), *(Q)}(0) = Q})(0) 
if and only if 

/ ([w£) x E((gojV)det^)E((PojV)det^Af)] ° J0)^ = 
Jm 

[ ([(cuD x E(Q)E(P)) o j'V; o j'0]detV>Mk 
Jm 
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but since this latter equation holds for all sections <j> of E it is equivalent to 

ojD x E((Q o j<i[>)deti> M )E((P o j^det^M) = [(ujD x E(Q)E(P)) o j^]d^ M 

up to a divergence. Using Lemma 13 1 U I II this is equivalent to 
uD x (dety M ^ E(Q) o ^)det^M^f (E(P) o^) = 

(wo jV)P*(E(Q)) o^)(E(P) o#)det^M 

or 

w[P x (detV>M^f)(E(Q) o jYO + detVM^f(^(E(Q)) o j^)D x ^ M ] 
(E(P) o#) = (wo^)(P>,(E(Q)) o #)(E(P) oj^) 

up to a divergence and where we have used equation (j2.1j) . (Notice that 
det^M cancelled from both sides of the equation since it is nonzero.) Finally 

dip 

since the last equation is true for all P and Q we must have DJdetipM—^ — ) = 

au 

dip 

and uo jtp = u(detipM) 2 (—R — ) 2 - (One may factor out (D X (E(Q)) a jip) over 

au 

appropriate domains from the left-hand side of the equation to verify that we 

dip 

must have DJdetibM—^ — ) = since Q is arbitrary, otherwise the equation 

ou 

may not hold for all P and Q. Also notice that by the arbitrariness of P 
and Q the divergence term must be zero.) Observe that the first condition 

dip 

implies that detipM—^ — is a constant, while the second condition was sim- 

ou 

plified using D x ipM = detipM since M is one-dimensional and so are its fibers. 



Theorem 3.1.1 Let E — ► M be a vector bundle with dim(M)=l, dim(E)=2 
and let the Poisson bracket {■,■} be defined by T> = ujD x . Then the in- 
duced transformation on the space of local functionals is canonical, i.e., 
{V(V), = V({V, Q}) for all V,Q E J 7 if and only if 

(i) D x (detip M ^^-) = 0, and 
(li) uojip = u;(det?/>M) 2 (^-^) 2 - 
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3.2 Case II: dim(M)=rc, dim(£)=m + n; V = u abl D l 

Suppose that M is n-dimensional with fibers of dimension m, and consider 
first order differential operators of the form T> = uj abl Di. Preserving the 
Poisson bracket: {V o^Qo V'}(0) = ({P, Q} o i/;)((f>) is equivalent to 

([uj abi DiE b ((Q o j^)det^/)E a ((P o j^)det^ M )\ o j(j>)u = 

M 

r ([(^AE fe (Q)E a (P)) o M o jfldetMv, 

M 

but since this latter equation holds for all sections <fi of E it is equivalent to 
o; a6i AE 6 ((go^)det^)E a ((PojV)det^) = [(^ M AE fe (Q)E a (P))ojV;]det^M 



up to a divergence. Using Lemma \3 . ( ) . 1 1 1 his is equivalent to 
c^A(det^M^f E d (Q) o#)det^M^|(E c (P) o j^) 



(u abi o jip)(Di(E b (Q)) ojip)(E a (P) o^)det^ 



M 



or 

dib d dib a 
^[A(det^ M ^f)(E d (g) o#) +detVM-^(D i E d (Q) oj^)D^ Mi 



()rs 'E c (P)ojip) = (^o^)(A(E b (Q))o^)(E a (P)o^) 



up to a divergence and where we have used equation (|2.1jl . Finally this can 
be rewritten as 

^'[^•(det^^f )(E 6 (Q) o jV) + det^M^f (A(E 6 (Q)) o #)^m] 



<9w c 



E a (P)o^) = (a; a&l o^)(A(E 6 (Q))o^)(E a (P)ojV) 



up to a divergence. As before for this to hold for all P and Q we must 

have Pjfdet^M = for all d,j satisfying uo cd i ^ for some c, and 

ou d 

, ,abi n • /, _ J» + ./, , ,cdj n ,/,« 
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Theorem 3.2.1 Let E — > M be a vector bundle with dim(M)=n, dim(E)= 
n+m, and let the Poisson bracket {■,■} be defined by V = u abl Di. Then the 
induced transformation \I/ on the space of local Junctionals is canonical, i.e., 
{V(V), *(Q)} = *({P, Q}) Jor all V,Q E J 7 ij and only ij 

dip b 

(i) uj cd ^ 7^ Jor some c =^> DAdetipM 5 ) = f or °^ ^ anc ^ 

air 

W u M oM = det^ M u cd ^^D^ M . 

3.3 Case III: dim(M)=l, dim(£)=2; £> = u/D/ 

Lastly we consider the case where we have a one- dimensional manifold M 
with one-dimensional fibers and n-th order differential operators of the form 
V = uj 1 Di = oj n D n + u) n ~ l D n ^i + • ■ ■ + a/-Di + u°. For simplicity we assume 
that / is a number rather than a multi-index since in this case the base 
manifold is one-dimensional so, for example, D n means (D x ) n and Dj = 
(Dx) 1 . Following similar analysis as before and skipping a few steps the 
Poisson bracket is preserved if and only if 

[u°(detrl> M ?p-E(Q) o jijj) + u I D I (deti> M ^E(Q) o j^)]de^ M 
ou Ou 

^(E(P) o#) = [(co°E(Q) + co I D I E(Q))E(P)] oj^det^M 
ou 

or 

[u°(de^ M ^E(Q)oj^) + f Q u I [D k {de^ M ^-)D I ^{E{Q)oj^)\ 
;E(P) oj^) = [(u; E(Q) + u/£>/E(Q))E(P)] o 
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up to a divergence. Now observe that using equation (J2.1)) we have 

Di(F o jip) = (DjFoj^ide^M) 1 + (D I _ 1 Fo J ^)(I,I-l) + 

(Dj_ 2 F o j^)(J, / - 2) + • • • + (DxF o j^)Dj^(det^ M ) 

where Dj = (D X Y and we use for the ( \ J ^ permutations of j 

detipM's and I — j D x 's with det"0M in the rightmost slot (this may be de- 



rived by induction). For example (3, 2) = D x ((detipM) ) + det'0M^a;(det , 
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This yields the conditions we must have for a canonical transformation: 

(^ oj^ = det^ (^) 2 . 

j 

{nW o^~tE({) o/Dj-^det^,^) (j, /); / = 1, 2, • • • , n. 

M j=j j=I V J / M 

Notice that we can combine condition (i) with condition (ii) if we let 
(0,0) = 1 and (j, 0) = for j ^ 0. To summarize 

Theorem 3.3.1 Let E — > M &e a vector bundle with dim(M)—l, dim(E)— 
2, and let the Poisson bracket {■,■} be defined by the n-th order differential 
operator!) = uj 1 Dj. Then the induced transformation ^ on the space of local 
functionals is canonical, i.e., {&(V),V(Q)} = ^({V, Q}) for all V, Q E T 
if and only if 

J ° i* = i£ E E ( j ) ^Dj-A'^m^M /); / = o, i, • • ■ , n. 

J=I j=I ^ ' 

dipE 

Remark If detip M and — — are constant then the above condition reduces 

ou 

to: J o ji, = <y (det^M) m (^f ) 2 ; I = 0, 1, • • • , n. 

4 Transformation of evolutionary systems and 
their conservation laws 

Recall that an evolutionary system of equations takes the form 

ut = V5H 

where D is a Hamiltonian differential operator and 7i the Hamiltonian of 
the system. The variational derivative S of a local functional V — J M Pv 
has the same components as E(P), the Euler-Lagrange operator applied to 
P (we use E(P) rather than E(Pz/) for simplicity). In the event of a canoni- 
cal transformation we have {\I/(P), ^f(Q)} = *&({V, <2}). Therefore if V is a 
conservation law for the evolutionary system u t = V5H with Hamiltonian H, 
then ty(P) is a conservation law for the evolutionary system u t = V5^(H) 
with Hamiltonian ty(H). 
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This transforms the evolutionary system and its conservation laws. We 
illustrate with a few examples. The first example is somewhat trivial but we 
find it a good point to start as an illustration of our approach. Also recall 
that the conditions for a canonical transformation are rather restrictive. In 
the second example, even though the action is orientation-reversing our con- 
clusions will still hold. The same applies for Example 3. 

Example 1 Consider the Korteweg-de Vries equation ut = u xxx + uu x with 

x 

the transformation if) defined by u = ku,x — — ,k > 0. Even though this 

/c 

transformation is somewhat trivial we find it useful to consider in this ex- 
ample. The Hamiltonian and Hamiltonian differential operator associated 

f 1 1 

with the KdV equation can be chosen to be H = / (—-u 2 + -u 3 )dx and 

dip 

V = D x respectively. Now observe that detipM—R — = 1 is a constant, and 

ou 

dip 

u o jip = uj(detipM) 2 (—^. — ) 2 = 1- So ip defines a canonical transformation 

of the Poisson bracket defined by T> = D x . This transformation transforms 

f 1 1 

our Hamiltonian to Ji! = / ( — -ui H — —u 3 )dx and consequently the KdV 

J 2k 4 x 6k 3 



equation to u t = -pru xxx + —uu x . 



1 1 

^3 Uxxx + p ' 



Notice that the conservation laws V\ — f -u 2 dx, V2 = f xu + -tu 2 dx, 

Jm 2 Jm 2 

and M. — / udx are trivially transformed to ^{V\) = / -u 2 dx^(V2) = 
Jm Jm 2 

j (— + -tu 2 )dx, and V(M) = / udx. 
Jm k 2 J M 

Example 2 Consider the Korteweg-de Vries equation 

Ui U xxx "T uu x 

as above but this time with the transformation defined by u = x 2 u,x = —. 

x 

Let M = R + and E = M x R 2 . In this case u x = (x 2 u x + 2xu)(— x 2 ) and 
the new Hamiltonian is given by H! = ty{H) = J (^x 6 u 2 x + 2x 5 uu x + 2x 4 u 2 — 

-x 4 u 3 )dx. So the KdV equation is transformed to u t = -D x E(iJ / ) where 
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H' = \x®u 2 + 2x 5 uu x + 2x 4 u 2 — |x 4 m 3 , i.e., we have 

u t = —2Ax 3 u — 2x 3 u 2 — x A uu x — 36x A u x — 12x 5 u xx — :r 6 ti xxx . 

The conservation laws V\ — [ -u 2 dx, V2 = [ xu + -tu 2 dx, and M. 



M 2 JM 



I udx are transformed to ^{V\) = f —\ ) x 2 u 2 dx^{V2)= [ ——+^-tu 2 dx, 
JM Jm 2 J M X 2 

*(M) = [ 
Jm 



and ty(M) = / udx. 

' M 



Example 3 The Boussinesq equation can be written as a system 

u t = v x , 
v t = \uu x + \u 



XXX • 



f 1 4 1 
It has Hamiltonian H = ( — u 2 H — u 3 H — v 2 )dx with associated Hamil- 

Jm 6 x 9 2 ; 
tonian differential operator 



V 





At 



Let M = R + and E = M x R 2 . Now consider the transformation u = 

x 2 u,v = x 2 v,x = — . The Hamiltonian H is transformed to H' = ^(H) = 

x 

f (A(-2x 3 u-x A u x ) 2 + -x 6 u 3 + -x A v 2 )^- = [ (--(Ax A u 2 + Ax 5 uu x + 

4 1 

a; 6 i4) + -x 4 m 3 + -x 2 v 2 )dx. This transforms the system to 
y ^ 

■u t = — 2xv — x 2 v x , 

1J-£ — $X U g 3y ^>xx ~3~*^' ^ 3*^ ^xxx' 

The Boussinesq equation has the conservation law 

p = /^*4 - 2mm ' + ^ 4 + 2uy2 - \ v D dx - 

This is transformed to = [ (6x 6 u 2 + 6x 8 u 2 x + -x 10 u 2 xx + \2x' uu x + 

Jm 6 

g 

2x 8 uu xx + 2x 9 u x u xx — 8x 6 u 3 + 8x 7 u 2 u x + 2x 8 uu 2 x + -x 6 u A + 8x 6 uv 2 + 8x 7 uvv x + 

1 

2x % uv 2 — 2x A v 2 — 2x 5 vv x — -x Q v 2 )dx. 
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